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1 Introduction 



It is known that flavor changing neutral currents(FCNC) induced rare B decays are very sen- 
sitive probes of new physics. The GIM suppression of FCNC amplitude is absent in many 
new physics scenarios beyond the Standard Model(SM), which could give large enhancement 
of FCNC processes over the SM predictions. To search such kind of signals is one of the 
most important goals of B projects BaBar, Belle, BTeV, and LHC-B. On the other hand, B 
rare decays also serve as laboratory for hadronic dynamics. Due to our poor knowledge for 
non-perturbative QCD, predictions for many interesting decays are always polluted by uncom- 
fortable large uncertainties, which have hindered us very much in extracting weak interaction 
information precisely from the available measurements. The well known example is two-body 
charmless B decays. It would be of great interesting to explore B rare decays which are induced 
by FCNC currents as well as involve few hadronic parameters. 

To the end, we will study the pure penguin annihilation decay — > (j) r y. Experimentally, 
Bd ~~ * ^7 ^ s ver y easy to be identified through the decay chain B® — > 07 — > (K~ K + )^, i.e, two 
charge tracks and one energetic photon, and the detecting efficiency should be high. However, 
if the decay rate is too rare, it would be very difficult to pick out the signals of the decay 
from its background, which comes from continuum (e + e~ — > qq with q=u,d,s) events with high 
energy photons originating from initial state radiation or e + e~ — > (7r rj)X with 7r r] — > 77pQ- 
CLEO Collaboration had performed a search for the decay, but found no evidence and put 
the upper limit B(B® — > 07 < 0.33 x 1CT 5 at 90% CL. To our best knowledge, there is no 
realistic theoretical study of B^ — > 07 in the literature. In this paper, we will investigate 
this decay. In naive factorization approach, this decay involve the simple matrix (0|s7 M s|O) 
and the same hadronic matrix element ( / -f\q / ~f^(l — / -f 5 )b\B q ) as the radiative leptonic decay 
which has been studied in Ref . [21 El HI E] with different framework. Beyond naive factorization, 
non-factorizable contribution should be included. The QCD factorization framework [H], which 
has been developed recently by Beneke et al., is employed to calculate 0(a s ) nonfactorization 
contributions arising from exchanging hard gluon between the two color octet currents s Q 7 M (l± 
7 5 )s^ and d^f(l - "f 5 )b a . We find B(B° d -> 7 ) = 3.6 x 10~ 12 in the SM. The decay rate is 
too small to be observed at the running B factories BaBar and BELLE. Any measurement 
of the decay at BaBar and BELLE would be the evidence of activity of New Physics. As an 
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example, we treat it as a probe of R parity violating couplings(RPV). Within the available 
lowest upper bound of RPV couplings, it is found that the branching ratio of the decay could 
be enhanced to 10~ 9 ~ 10~ 8 , which might be measured at LHC-B, BTeV and the planning 
super high luminosity B factories at KEK and SLAC. 



The — > 07 decay is illustrated schematically in fig.l. which is dominated by the photon 
radiating from the light quark in the B meson. Obviously, the amplitude of fig.l is suppressed 
by power of AQco/n^b because <fi meson must be transversely polarized and B meson is heavy. 
The diagrams with the photon radiating from the heavy b quark and the energetic strange 
quarks of the <j) meson are suppressed by power of A.Q CD /ml, which will be neglected in this 
paper. The situation is similar to that of the annihilation contributions in B — > K*^ decaysjZj. 



Figure 1: The leading diagram for B® — > cjyy. Radiation of the photon from the remaining three quark 
lines are further suppressed by power of A^c^/mfe. The symbol ® denotes the insertion of penguin 
operators O3_io- 

We begin our study with the effective Hamiltonian in the SM relevant to the decay [S] 




07 in the Standard Model 





(i) 



For convenience, we list below the operators in Tt e ff for b — > dss: 
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daj^Lba ■ s^^Rsp , 6 
-d^Lba ■ esSpipRsp , 0$ 



d a YLb p 



5 



7 



\d a ^Lbp ■ esSpJuRs 



(2) 



a 



9 



10 
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Where a and (3 are the 577(3) color indices and L — (1 — 7s)/2, i? = (1 + 7s)/2. The Wilson 
coefficients evaluated at /i = m?, scale arejH] 



C a = 0.014, 



-0.035, C 5 = 0.009, 



-0.041, 



(3) 



C 7 = -0.002/137, C 8 = 0.054/137, C g = -1.292/137, C 10 = 0.263/137, 
Using the effective Hamiltonian and naive factorization hypothesis, it is easy to write down the 



amplitude for B d 



Gp 
'V2 



v tb v; d 



9i 

Nr. 



N 



N 



x v^^m^{-W^± elW + i Kvl ■ <L){v ■ q) - (t£ ■?)(«• e±)] j(4) 

where 7/^_ and are transverse polarization vectors of and photon respectively. The form 
factor F v is defined by0 El H El 



( 7 (e*, g)| J 7m (1 - 75)&|5°) = v 7 ^: f-Fye^e'W + iF A (e>i; - g^-e*) 



(5) 



To leading power of 0(1/ mb), Fy and Fa readj21El 

Fv 



2V2E^ 



(6) 



The contributions of strong penguin operators arising from the renormalization group evolu- 
tion from the scale \x = Mw to \i = nib is very small due to the cancellations between them: 
C3 ~ — C4/3 and C5 ~ — Cq/3. Obviously the amplitude is dominated by electro- weak pen- 
guin(EWP). 

Now we can write down the helicity amplitude 



M+, 



G 

I -J= V th V* d \JA~Kot e F vf^nifyM-B 



«3 + a 5 - -{a 7 + a 9/ 



0. 



(7) 
(8) 



where a, = Ci + Ci + i/Nc- It is interesting to note that the <fi meson and the photon in the 
decay are right-handed polarized in the SM. It is also easy to realize that the decay is very 
rare because of helicity suppression as well as small V t d, fs and f^. Using = 254MeV|H|, 
f B = 180MeU, \Vtd\ = 0.008, N c = 3 and A^ 1 = / dl + <$> B1 = (0.35GCU)- 1 Ej, we get 



B(B° d -> 7 ) = 3.5 x 10" 13 . 



(9) 



In the above calculations, non-factorizable contributions are neglected. However, the leading 
non-factorizable diagrams in Fig.2 should be taken into account. For this purpose, the QCD 
factorization framework[6j invented recently by Beneke, Buchalla, Neubert and Sachrajda is 
very suitable. The framework incorporate many important theoretical aspects of QCD like color 
transparency, heavy quark limit and hard-scattering which allows us to calculate nonfactorizable 
contributions systematically. 

To calculate the nonfactorizable diagrams as depicted by Fig.2, we take photon and <fi 
meson flying along n_ = (1, 0, 0, —1) and n + = (1, 0, 0, 1) directions respectively. We need the 
two-particle light-cone projector for B meson and meson 



M °e = w/ bMb { {1+ ths[*Bi(i+)+ rl-$B2(! + )]} a , (io) 

Mi h = -^{^? ) ( M )+, WCTe ix<7^5^} 5 , (ii) 



which encode the relevant non-perturbative bound state dynamics of the initial B meson and 
the final <p meson. <&bi{1+) and §32(1+) are the leading twist light-cone distribution function 
of B meson [TU]. g±\u) and g±\u) are twist-3 distribution amplitude of <p meson[5], <7j* I 
dg^f' \u) I 'du. The detail discussions on these projectors could be found in Refs. |§l ITU1 ITT] . 
The amplitude for B® — > 07 can be presented by the schematic formula 



u 



(<h\He ff \B) = 



[a' 3 + a' 5 ) - - (a' 7 + a' 9 ] 



x v/tel/^m^j-e^^elW + i [(vl ■ e* ± vq) - (vl • q)(v • e*)]}{12) 
The 0(a s ) corrections are summarized in a'^s which are calculated to be 

a 3 = a 3 + - — j-C 4Fl , (13) 

a 5 = a 5 + -—j-C 6 F 2 , (14) 

, a s C F f<f> n „ H K x 

a 7 = a 7 + - — -C B F 2 , (15) 

a 9 = a 9 + - — -C l0 F u (16) 

where F\ 2 arise from one gluon exchange between the two currents of color-octet penguin 
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Figure 2: Nonfactorizable diagrams for — ► ^7. Other diagrams with the photon radiating from 
the heavy b quark and the energetic quark lines are suppressed. 



operator 04,6,8, 10 as shown by fig. 2 
Fi = 



\ 4 



-14-3i7T- 12 In — 



u 



lnw 



u 



7T 

y 



2Li 2 { 



u 



u 



(17) 



1 . I («)/ \ . 



-14 - 3i?r - 12 In 



+ 5 



— In it 1 

it / 3 



2Lz, 



,u- 1. 



(18) 



In the calculation, the MS renormalization scheme is used. We have neglected the small effect of 
box diagrams and the diagrams with photon radiating from energetic strange quarks, which are 
further suppressed by Aqcd/M b . We also have neglect l\ terms entered in the loop calculation 
which are the higher twist effect, in this way, the integral involved 3>b2{1+) absents and the 
remain integrals are related to the form factor Fy. 

Including 0(a s ) contributions, the helicity amplitude is 

.Gf ^ 



At 



1 -j= Vtb y t *d \J^a e F v f^m^MB 



a', + ai 



= 0. (19) 

To give numerical result, we take \x = mj and = 215MeV^. The branching ratio is 
estimated to be 

B{B° d -> 7 ) = 3.6 x 1(T 12 . (20) 

It indicates that the decay is too rare to be measured at the running B factories BaBar and 
BELLE. However, the decay may be accessible at LHCB and the planning super high luminosity 
B factories at KEK and SLAC. Furthermore, it could be enhanced by New Physics, and the 
enhancement might be large enough for being measured at these facilities. 



3 —> 07 as a probe of R parity violation 

In the minimal supersymmetric standard model( MSSM ) ^2] a discrete symmetry called R- 
parity is invoked to forbid gauge invariant lepton and baryon number violating operators. The 
R-parity of a particle is given by [12] R p = (— i) L + 25 + 3 - B ) where L and B are lepton and baryon 
numbers, and S is the spin. However there is no deep theoretical motivation for imposing 
R-parity and it is interesting to explore the phenomenology of R-parity violation [I3|. 
We start our exploration from the R-parity violating superpotential 

W$ = -XLLE C + \'LQD C + -\"U C D C D C . (21) 

We are interested in the A' and A" terms since they are relevant to the process b — > dss. Writing 
down all indices explicitly, it reads 

(22) 

where U? a , and D k ^ are the superfields of the right handed quarks and/or squarks respec- 
tively, the superscript c denotes charge conjugation, i, j and k are generation indices and a, 
(3 and 7 are SU (3) color triplet indices. It follows from the antisymmetry of e al31 that A 'Lu is 
antisymmetric in the last two indices. 

In terms of four components Dirac spinor, from W% we can read 

+ y ijk v Li d k Ldj + h.c. (23) 
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Wv = e^X' ijk L ia Q jba D c k( , + -e a ^ 



\" TT C D c n c 
A i\ik\ {J ici U jp u k" 1 



Kjk u ma (d c jl3 Rd kl - {j <-+ k} 



From C e ff, we get the effective Hamiltonian for b —>■ dss 



n 



ml 



. e <xPl e a P'l' 



X: 23 X^(s c M(s^Ld%) + X^X^Rb^dyLsV, 



^32^/21(^^7) (dyLs c g,] 



in? 



KaiKzi&aRsctXdpIAp) + X'* 13 X' i22 (s a Ls a )(dpRbp) 
+Kz2Ki2(daRs a ){s p Lbp) + X'* 2Z X' i21 {d a Ls a ){spRbp) 

Contracting e^e 01 ^' 1 ' and performing Fierz transformations, we get 

2 



H$ (b — > dss) 



m 2 "m~ii2 



> (3 7 AtJ Rs /3 )(ci 7 7 M J R6 7 ) - (s^^Rs^^d^ Rb^ 



2ml 



(24) 



KzxKA^l^h^dp^Rso) + X'* lz X[ 22 {s alll Rbp){d^Ls c 



-KwKiiidctl^Lb^is^Rso) + Af 23 X 21 (d Q 7 /i i?6 /3 )(s /3 7^Ls a ) 



(25) 



where we have used the relation 



{q'^Rq) . 



(26) 



The effective Hamiltonian is near to the idea form. Using renormalization group to run it from 
sfermion mass scale (lOOGeV assumed ) down to nit, scale, it reads 

2 

TC% (b — > dss) = 



2 '/ "i23"il2 
" L Ui 



V^mKu {spl,Rsp){d^Rb,) - {sfnrRsJfa-fflbt 



2ml 



V 



\nK22{ s al^Lbp){dp^Rs a ) + X* n X i22 {s a -i^Rbp){d^Ls c 



+X' i32 X'* 2 (d a -f^Lb l3 )(sf 3 ^Rs a ) + X% z X' i2l {d all ,Rbp){s^Ls a ) 

™ th v = and P° = 11 ~ 3 n /- 

Now we are ready to write down R-parity violation contributions in — > 07 decays 



(27) 



M# (B° d - 7 ) = V^eFvUm^ [i{ v l ■ el) (v q) + e^fe^q 



1 A' A'* «- 8 /* ( 1 /^^"ll 
8? V ^ \n~ c + ^W c F2 )) 



y/A7ra' e F v f (j) m (f) ]i(r]* ± ■ e* ± )(v ■ q) - e l , ua(3 7]*/e*^v a q 13 
x — A' A'* n-W° (— + ——F,\ 



(28) 



S 



We note that only u channel squarks and sneutrino mediated terms contribute because of 
(7|d(l ± 75) &| Bj) = 0. Immediately, we derive 




RPV Lambdas 



Figure 3: The branching ratio of — > (fry as a function of R-parity violation couplings lA-'^A-'j^Kthe 
thick solid curve), |A- 32 A-i 2 | (the thin solid curve) and | X' i2 i A" 2 * 3 | (the dash curve) respectively. The 
horizontal line is the SM prediction. 

To present numerical results, we assume that only one sfermion contribute one time and 
universal mass lOOGeV for sfermions um. Our results for R parity violation contributions are 
summarized in fig. 3. We note that the decay is very sensitive to R parity violating couplings. 
Within the available lowest upper bounds for |A" 23 A" 21 | [TBI FT01 fTTj . 

IA^ASlI < 6 x ID" 5 (^) 2 , |A:3 2 A;* 12 | < 4 x ID" 4 (^) 2 , (31) 

B(B® — * 07) could be enhanced to 10 -9 ~ 1CT 8 which may be inaccessible at BELLE and 
BaBar. However, it is large enough for LHCb, BTeV and the planning super high luminosity 
B factories at KEK and SLAC. 
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4 Conclusion 



We have studied pure penguin radiative annihilation process B® — > 07 by using QCD factor- 
ization for the hadronic dynamics. We find that non-factorizable contributions are larger than 
factorizable contributions in — > 07 decays. We estimate B{B^ — > 07) = 3.6 x 10~ 12 in 
the SM. The smallness of these decays in the SM makes it a sensitive probe of flavor physics 
beyond the SM. To explore the potential, we have estimated the contribution of RPV couplings 
to this decay. We have found that the decay B® — > is very sensitive to the X"U C D C D C and 
X'LQD terms in the R-parity violating superpotential. Within the upper limits for A'/^Af^ and 
A- 32 A^ 2 , it is found that B(B% — > 07) could be enhanced to order of 10~ 9 ~ 10~ 8 . We also 
note that the meson and the photon are right-handed polarized in the SM, but they can be 
left-handed polarized in RPV supersymmetry. We find that A£ 32 A^ 2 RPV couplings contribute 
to the right-handed polarized magnitude while A-^A-'j^ and A£ 21 A£ 23 RPV couplings responsible 
to the left-handed polarized magnitude. Theoretically, measuring the polarization of photon 
(or the meson ) can be used to probe the relative strength of RPV couplings. However, the 
interesting phenomena is faded by the small branching ratio of the exotic decay. In literature, 
there are interesting proposals for probing new physics by measuring the photon polarization 
in radiative decays B — > if*7|18]. It always needs large amount of experimental data. For the 
more exotic decay B^ — > 07, this phenomena will be remained academic unless the continuum 
background could be suppressed efficiently in data analyzing at the future B facilities. The 
interesting decay may be too rare to be accessible at Babar and Belle. However, the decay 
rate might be studied by LHCB at CERN, BTeV at Fermilab and the planning super high 
luminosity B factories at KEK and SLAC to probe new physics. 
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